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Abstract The hadronic kj^ -spectrum inside a high energy jet is determined including corrections of relative magni- 
tude O (-yol) with respect to the Modified Leading Logarithmic Approximation (MLLA), in the limiting spectrum 
approximation (assuming an infrared cut-off Qo — A ) and beyond (Qo 7^ A ). The results in the limiting 

spectrum approximation are found to be, after normalization, in impressive agreement with measurements by the CDF 
collaboration. A new integral representation for the "hump-backed" plateau is also reported. 

PACS. 12.38.-t Quantum chromodynamics - 12.38.Bx Perturbative calculations 
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The production of jets - a collimated bunch of hadrons - in 
e+e~, e^p and hadronic collisions is an ideal playground to 
investigate the parton evolution process in perturbative QCD 
(pQCD). One of the great successes of pQCD is the existence 
of the hump-backed shape of inclusive spectra, predicted in [ 1] 
within MLLA, and later discovered experimentally (for review, 
see e.g. [2]). Refining the comparison of pQCD calculations 
with jet data taken at LEP, Tevatron and LHC will ultimately 
allow for a crucial test of the Local Parton Hadron Duality 
(LPHD) hypothesis 1 3 1 and for a better understanding of color 
neutralization processes. 

The inclusive fc^ -distribution of particles inside a jet has 
. been computed at MLLA accuracy in the limiting spectrum ap- 
' proximation f?], i.e. assuming an infrared cutoff Qo equal to 
"^QCD = ^^(^o/^qcd) " ^ review, see also [5|). 

MLLA corrections, of relative magnitude O (\/a7) with re- 
spect to the leading double logarithmic approximation (DLA), 
were shown to be quite substantial [4|. Therefore, it appears le- 
gitimate to wonder whether corrections of order O (as), that is 
next-to-next-to-leading or Next-to-MLLA (NMLLA), are neg- 
ligible or not. 

The starting point of this analysis is the MLLA Master 
Equation for the Generating Functional (GF) Z = Z{u) of 
QCD jets [5 1, where u = u{k) is a certain probing function 
and fc, the four-momentum of the outgoing parton. Together 
with the initial condition at threshold, the GF determines the 
jet properties at all energies. For instance, the single inclusive 
spectrum can be derived from the GF by differentiating with 
respect to u = u{k), and the solution of the equations can be 
written as a perturbative expansion in -Ja^. At high energies 



this expansion can be resummed and the leading contribution 
be represented as an exponential of the anomalous dimension 
7(q!s). Since further details to this logic can be found in L5][6l, 
we only give the symbolic structure of the equation for the GF 
and its solution, which we write respectively as 



ay 



exp 



where ^{ag) can be expressed as a power series of -Ja^ 



^{a,{v'))dy'\, (1) 



(2) 
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In this logic, the leading (DLA, 0{y/oi^)) and next-to-leading 
(MLLA, 0{as)) approximations are complete. The next terms 

(NMLLA, 0(0:^^^)) are not complete but they include an im- 
portant contribution which takes into account energy conserva- 
tion and an improved behavior near threshold. Indeed, some re- 
sults for such NMLLA terms have been studied previously for 
global observables and have been found to better account for 
recoil effects and to drastically affect multi-particle production 

nam. 

The equation in ^ applies to each vertex of the cascade 
and its solution represents the fact that successive and indepen- 
dent partonic splittings inside the shower, exponentiate with 
respect to the evolution-time parameter dy = dO/O, where 
6* <C 1 is the angle between outgoing couples of partons. The 
choice of y follows from Angular Ordering (AO) in intra-jet 
cascades. It is indeed the suited variable for describing time- 
like evolution in jets. Thus, Eq. ([T]) incorporates the Markov 
chains of sequential angular ordered partonic decays which are 
singular in and 7(0:5) determines the real rate of inclusive 
quantities growth with energy. 

While DLA treats the emission of both particles as inde- 
pendent by keeping track of the first term ^ \fa.s ™ without 
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constraint, the exact solution of the MLLA evolution equation 
(partially) fulfills the energy conservation in each individual 
splitting process (z + {1 — z) — 1) by incorporating higher 

order (as ,n>l) terms to the anomalous dimension. Sym- 
bolically, the first two analytical steps towards a better account 
of these corrections in the MLLA, NMLLA evolution, which 
we further discuss in |[9l, can be represented in the form 



{as + as£ ^ In z)dz ~ + 



^3/2 



(3) 



where £ = ln(l/a::) ~ as with a; ^ 1 (fraction of the 
jet energy taken away by one hadron), z ^ 1 for hard partons 
splittings such as g ^ qq. . . (this is in fact the region where the 
two partons are strongly correlated). 

Energy conservation is particularly important for energetic 
particles as the remaining phase space is then very limited. On 
the other hand, a soft particle can be emitted with little impact 
on energy conservation. Some consequences of this behavior 
have also been noted in 1 10|: 

(i) the soft particles follow the features expected from DLA; 
(a) there is no energy dependence of the soft spectrum; 
(Hi) the ratio of soft particles r = Ng/Nq in gluon and quark 
jets is consistent with the DLA prediction Nc/Cp = 9/4 
(see the measurement by DELPHI lITTl '). 

The present study makes use of this logic (|3]l to evaluate 
NMLLA contributions to the single inclusive fc^ -distribution. 
The main results of this work have been published in |9 12 |. 
Experimentally, the CDF collaboration at the Tevatron reported 
on fc^ -distributions of unidentified charged hadrons in jets pro- 
duced mpp collisions at = 1.96 TeV Iil3il . 



2 MLLA evolution equations 



We start by writing the MLLA evolution equations for the frag- 
mentation function {x / z, zEOq^Qq) of a parton B (en- 
ergy zE and transverse momentum = zEOq) into a gluon 
(represented by a hadron h (energy xE) according to LPHD 
|3J) inside a jet A of energy E for the process depicted in Fig.[il 
As a consequence of angular ordering in parton cascading, par- 




A = (Q, G) 



Figure 1. Parton A with energy E sphts into parton B (respectively 
C) with energy zE (respectively, (1 — z)E) which fragments into a 
hadron h with energy xE. 



tonic distributions inside a quark and gluon jet, Q, G(2;) = 
x/ z Dq q (x/z, zEOq, Qo), obey the system of two coupled 



equations ||6l (the subscript y denotes d/dy) following from 
the MLLA master equation described by (|T) 



Qy 



Gy — 



dz^mz) 



l-z)-Q]+G{z) 



(4) 



az — 

TT 



<1>1{z){1-z)(g{z) + G{1-z)-G) 



G 



(5) 



^^(z) denotes the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) splitting functions 15], 

<P^q{z)^CF(^-^ + , f^(z)(l - z) = 2iV,Q + <t>l{z)^ . 

<Pliz)=Tn[z' + il-zr], 

where = z - 2 and = (z - 1)(2 - z(l - z)) are 

regular as z ^ 0, Gf = iN^~l)/2Nc, Tr = nf/2 (N^ = Sis 
the number of colors for SU{3)c and nf — Sis the number of 
light flavors we consider). The running coupling of QCD (a^) 
is given by 



27r 



ANcMi + y + X)' 



00 = ^ '-^N.-'-Tr 
^ ' 4Nc\ 3 3 



and 



e^{l/x) ,y = \n{kjQo) , X = In {Qo/ A^^J , 

where Qo is the collinear cut-off parameter Moreover, 

G=G{l) = xD'^{x,EOo,Qo), 
Q = Q{l)^xD^{x,EOo,Qo)- 

At small a; <C z, the fragmentation functions behave as 



B{z) w Pb ln-,ln— — 
x Qo 



p% (lnz + £,y)^ 



Pg being a slowly varying function of two logarithmic vari- 
ables ln(z/x) and y that describes the "hump-backed" plateau. 
Since recoil effects should be largest in hard parton splittings, 
the strategy followed in this work is to perform Taylor expan- 
sions (first advocated for in 1 151) of the non-singular parts of 
the integrands in ( I4l5l l in powers of Inz and ln(l — z), both 
considered small with respect to £ in the hard splitting region 
z~ 1-z = 0(1) 

B{z) = B(l) + B^(l)lnz + (In^z) ; z 1 - z . (6) 

Each ^-derivative giving an extra y^oj factor (see |l6l), the terms 
Bi{l)lnz and 5^(1) In (1 - z) yield NMLLA corrections to 
the solutions of (|5]). Making used of (|6]), after integrating over 
the regular parts of the DGLAP splitting functions, while keep- 
ing the singular terms unchanged, one gets after some algebra 
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(72 = 27V,a,/7r) miia 



N, 



G{£, y)^5{£) + I'm' Cdy'-flil' + y') \l - a^5{e - 1) 



+a2Sie' -e)M^',y') G{e',y'), 



(7) 



(8) 



with y) = Ge{£, y)/G{i, y). The MLLA coefficients 
ai =3/4 and ai « 0.935 are computed in fS) while at NM- 
LLA, we get 



7 Gf /5 7r2\ ^ _ 



(9) 



67 IZufTnGF 
-=36-T"l8^]€"°-'^'- ^^'^ 

Defining ?/) = 7o + y)G{£, y), we can exactly solve the 
self-contained equation ^ by performing the Mellin transform 



F{e,y)^ 



dojdv 



(11) 



(27ri)2 

Inserting (fTTl i into (O we obtain, after some algebra, the differ 
ential equation 



where the integration is performed with respect to t defined by 

1 y 

a = - m — + ir and with 



^B{T,y,l) 



Z{T,y,£) 



cosh a — — sinh a 

y + i 

t + y a 



Po sinh a 



B/2 



lB{2VZ{T,y,e)), 



e- 



y 



a 



f3o sinh a 



cosh a 



y + i 



■ sinh a 



B = ai/fio and Ib is the modified Bessel function of the 
first kind. To get a quantitative idea on the difference between 
MLLA and NMLLA gluon inclusive spectrum, the reader is re- 
ported to the appendix B of |9|. The magnitude of 02, however, 
indicates that the NMLLA coiTections to the inclusive quark 
jet spectrum may not be negligible and should be taken into 
account. After solving ([H]), the solution of (|7| reads 

y)=^ [g(^, y) + (ai - ai) G,{1, y) (15) 
+ {ai (ai - ai) + ^2 - 02) Gu{l, y) + 0(7^). 

It differs from the MLLA expression given in |4| by the term 
proportional to Gu, which can be deduced from the subtraction 
of (Cf/A^c)x® toEq. ^. 



/3o A.F- — - — =- + ai- - a^-T. (12) 

\ oll) av j V uv v V 

Finally, after inserting the solution of (fTSl i into (fTTT i. the new in- 
tegral representation for the inclusive spectrum (hump-backed 
plateau) with NMLLA accuracy reads 



G(^,2/)=(^ + y + A) 



ijj(y ^ s) 



dujdv 
(2^ 



{uj + s)v 



where 



/3o(^- 



ai 



'% (13) 



— {uj-iy), (73 = — + A. 
Po Po 



However, computing ( fT3] l numerically is quite a challenging 
task. As can be seen, the NMLLA coefficient 02 is very small. 
This may explain a posteriori why the MLLA "hump-backed 
plateau" agrees very well with experimental data 1 1 , 14|. There- 
fore, the NMLLA solution ( fTsT i of (O can be approximated by 
the MLLA solution of G{(, y) {i.e. taking 02 ~ 0), which 
will be used in the following to compute the inclusive k^- 
distribution. As demonstrated in |6|, taking the limits a2 ~ 
and A « in ( fT3] l. the integral representation can be reduced to 
the known MLLA formula Q 



G{i,y) = 2 



/3o 



2 dr 



-Ba 



^B{T,y,i), 



(14) 



3 Single inclusive ^_L-distribution of 
charged hadrons in NMLLA 

Computing the single inclusive k±— distribution requires the 
definition of the jet axis. The starting point of our approach 
consists in considering the correlation between two particles 
(hi) and (h2) of energies Ei and E2 which form a relative 
angle inside one jet of total opening angle Oq > lfT6l . 
Weighting over the energy E2 of particle (h2), this relation 
leads to the coiTelation between the particle (h=hl) and the 
energy flux, which we identify with the jet axis (see Fig. O 
||4j|. Thus, the correlation and the relative transverse momentum 
k± between (hi) and (h2) are replaced by the coiTelation, and 
transverse momentum of (hi) with respect to the direction of 
the energy flux. Finally, we obtain the double differential spec- 
trum d^N/dx d0 of a hadron produced with energy Ei — xE 
at angle (or k± « xE0) with respect to the jet axis. As 
demonstrated in |4 |, the correlation reads 



dx d In 6» d In 6> 



F'X, ix,0,E,0o) 



(16) 



where F^l^ is given by the convolution of two fragmentation 
functions 

Pao=J2 I duDi^{u,E0o,uE0)D'x(-,uE0,Qo 



A=g,q 



(17) 
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u being the energy fraction of the intermediate parton A. 
describes the probabiHty to emit A with energy uE off the par- 
ton Ao (which initiates the jet), taking into account the evolu- 
tion of the jet between 0q and O. D\ describes the probabiHty 
to produce the hadron h off A with energy fraction x/u and 
transverse momentum k, uEO > Qo (see Fig.|2|. As dis- 



{C)g^q can be seen as the average color charge carried by the 
parton A due to the DGLAP evolution from to A. Introduc- 
ing the first and second logarithmic derivatives of D\, 




(Jet Axis) 



Figure 2. Inclusive production of hadron h at angle inside a high 
energy jet of total opening angle 0o and energy E. 



D'Xii + lnu,y) ''w" + In u, Ye + In m) , 



(18) 



with Y0 = £ + y = IhEO/Qq. Performing the Taylor ex- 
pansion of D to the second order in (In u) and plugging it into 
Eq. ([T7]i leads to 



^Pao- J2 [ duuDi^iu,EOo,uE0)D'Xii,y) (19) 
+ f duu\nuDi^iu,Eeo,uEO)'^^'^^^'y^ 



A=g,q 



+- I duulT? uD^^{u,E0o,uE0) 

2 ^ J 



A=g,q 



d'D'X{l,y) 

dp ■ 



Indeed, since soft particles are less sensible to the energy bal- 
ance, the correlation ( [TtI i is proved to disappear for these par- 
ticles, leading to the sequence of factorized terms written in 
( fT9] l. The first two terms in Eq. (fT9] ) correspond to the MLLA 
distribution calculated in [4| when I?^ is evaluated at NLO 
and its derivative at LO. NMLLA corrections arise from their 
respective calculation at NNLO and NLO, and, mainly in prac- 
tice, from the third line, which was computed in |9 12J. In- 
deed, since x/u is small, the inclusive spectrum D^{£, y) with 
A = G,Q are given by the solutions ( fT3] l and ( fTSI l of the next- 
to-MLLA evolution equations (|7]i and dHJ respectively. How- 
ever, because of the smallness of the coefficient 02 (see (fTOli). 
G{£, y) shows no significant difference from MLLA to NM- 
LLA. As a consequence, we use the MLLA expression ( fT4l i 
for G{£, y), and the NMLLA ^ for Q{t, y). The functions 
and F^^ are related to the gluon distribution via the color 
cuiTents {C)g,q defined as: 



xF^ 



N, 



G{l,y). 



(20) 



i'>pA,i + '>pA,ie)i£,y) 



1 dP A{i,y) 

blil.y) dl 

1 d^DA{l,v) 
D\ie,y) dp 



= 0{as). 



Eq. ( fT9l ) can now be written as 



A=g,q 

+hu\n^ u)X,{ipA/ + '4'A,ee){i,y) 



with the notation 



cussed in |]4l, the convolution ( fTTI l is dominated by u ^ 1 and 
therefore D^^ (u,E0o,uE0) is determined by the DGLAP 
evolution jS). On the contrary, the distribution D'^ (^^,uE0, Qo) 
at low a; <C u reduces to the hump-backed plateau. 



>Ao 



[ du {u In' u) Di^ {u, E0O, uE0) 
Jo 

I du (u In* u) {u, E0O, E0) . 
Jo 



The scaling violation of the DGLAP fragmentation function 
neglected in the last approximation is a NMLLA coiTection to 
< u >. It however never exceeds 5% |9| of the leading con- 
tribution and is thus neglected in the following. Using ( |20] |. the 
MLLA and NMLLA corrections to the leading color cuiTent of 
the parton Aq = g-.q read 



+ Cf(u In 



S: / iNiv 
^{C^/Ao 



NMLLA-MLLA 



(21) 

= N,{u\ii'u)%^{i;l, + i;g,u) (22) 
+ Cf {uln^ u)\ [li^lj, + 



The MLLA correction, O (y^oj), was determined in |^ and 
the NMLLA contribution, O (a^), to the average color current 
is new. The latter can be obtained from the Mellin moments of 
the DGLAP fragmentation functions 

2?1(J,0= f' duu^-'DX{u,0, 
Jo 



leading to 



(uln^u)i^ 



dj 



^vX{j,aE0o)-aE0)) 



(23) 



i=2 



Plugging ( |23] ) into the resummed expression < C >g.q of the 
color current as written in ( l20l l. these quantities for gluon and 
quark jets are determined analytically |9|. For illustrative pur- 
poses, the LO, MLLA, and NMLLA average color current of a 
quark jet with Fg^ = 6.4 - coiTesponding roughly to Tevatron 
energies - is plotted in Fig.[3]as a function of y, at fixed £ = 2 
for A « 0. As discussed in [4], the MLLA corrections to the LO 
color current are found to be large and negative. As expected, 
the correction O {as) from MLLA to NMLLA proves much 



Redamy Perez Ramos: Transverse momentum spectra of charged hadrons in jets at the Tevatron 



5 



— LO 

- - - MLLA 
NMLLA 



- MLLA 

- NMLLA 



Figure 3. The color current of a gluon jet (left panel) and a quark jet 
(right panel) with Ye^ = 6.4 as a function of y at fixed £ = 2 for 
A « 0. 




Figure 4. NMLLA corrections to the color current of a quark jet with 
Yeg ~ 6.4 and ^ = 2 for various values of A. 



smaller; it is negative (positive) at small (large) y. This calcu- 
lation has also been extended beyond the limiting spectrum, 
A 7^ 0, to take into account hadronization effects in the pro- 
duction of "massive" hadrons, m — O [Qo) 1 18 1. We used, ac- 
cordingly, the more general MLLA solution of dSjl with 02 = 
for an arbitrary A 7^ 



G{l,y) = {£ + y + X) 



U!{v -\- S) 



^ 



V -\- S 



ly + s 



from which an analytic approximated expression was found us- 
ing the steepest descent method ifTTl ; <72 = and 0-3 = A have 
been set in ( fTST l. However, G{£, y) is here determined exactly 
from an equivalent representation in terms of a single Mellin 
transform (which reduces to ( fT4] l as A ^ 0) 1 18 1 



£ + y + X 



duj 



/3o B (B + 1) 2Tn 



(24) 



X <P{-A + B + l,B + 2, -Lo{l + y + \)) K,{uj, A) 

which is better suited for numerical studies. The function K. 
appearing in Eq. (l24l i reads 



/C(w,A) 



r{A) 
r{B) 



[lo A)^!f(A,B + l,cjA), 



(25) 



where A — l/(/3o lo), B — ai/f3o, and <P and tf' are the con- 
fluent hypergeometric function of the first and second kind, re- 
spectively. The NMLLA (normalized) corrections to the MLLA 
result are displayed in Fig.|4]for different values A = 0, 0.5, 1, 
by using (l24l i. It clearly indicates that the larger A, the smaller 
the NMLLA corrections. In particular, they can be as large as 
30% at the limiting spectrum (A = 0) but no more than 10% 
for A = 0.5. This is not surprising since A 7^ (Qq 7^ ^qcd-* 
reduces the parton emission in the infrared sector and, thus, 
higher-order corrections. The double differential spectrum 
(d'^N/di dy), Eq. (fTSI l. can now be determined from the NM- 
LLA color currents (|22] | using the MLLA quark and gluon dis- 
tributions. Integrating it over £ leads to the single inclusive y- 
distribution (or -distribution) of hadrons inside a quark or a 
gluon jet: 



diV 
dy 



9,1 



dN 
dk~ 



9,g 



d£ 



d^N 
d£dy 



(26) 



where 



d^N 



didyj, + 

(27) 

is the explicit expression of the double differential distribu- 
tion ( fTSI l. The MLLA framework does not specify down to 
which values of £ (up to which values of x) the double differ- 
ential spectrum {d'^N/d£dy) should be integrated over Since 
{d'^N/d£dy) becomes negative (non-physical) at small val- 
ues of £ (see e.g. [4J), we chose the lower bound ^min so as 
to guarantee the positiveness of {d'^N/d£dy) over the whole 



< £ < Yog range (in practice. 



1 and 



2). 



Having successfully computed the single kj^ -spectra including 
NMLLA corrections, we now compare the result with existing 
data. The CDF collaboration at the Tevatron recently reported 
on preliminary measurements over a wide range of jet hard- 
ness, Q = E0Q, in pp collisions at ^/s = 1.96 TeV |13|. 
CDF data, including systematic errors, are plotted in Fig. [Alto- 
gether with the MLLA predictions of |U| and the present NM- 
LLA calculations, both at the limiting spectrum (A = 0) and 
taking A^^^ — 250 MeV; the experimental distributions suf- 
fering from large normalization errors, data and theory are nor- 
malized to the same bin, ln(fc^/lGeV) — —0.1. The agree- 
ment between the CDF results and the NMLLA distributions 
over the whole fc^^ -range is particularly good. In contrast, the 
MLLA predictions prove reliable in a much smaller k^^ inter- 
val. At fixed jet hardness (and thus Ysg), NMLLA calcula- 
tions prove accordingly trustable in a much larger x interval. 
Despite this encouraging agreement with data, the present cal- 
culation still suffers from certain theoretical uncertainties, dis- 
cussed in detail in |9|. Among them, the variation of A 

' QCD 

-giving NMLLA corrections- from the default value A = 

° QCD 

250 MeV to 150 MeV and 400 MeV affects the normahzed k^- 
distributions by roughly 20% in the largest ln(fcj^/l GeV) = 
3 GeV-bin at Q = 100 GeV. Also, cutting the integral (|26ll at 
small values of £ is somewhat arbitrary. However, we checked 
that changing £min from 1 to 1.5 modifies the NMLLA spec- 
tra at large k^ by ~ 20% only Q. Finally, the kt-distribution 
is determined with respect to the jet energy flow (which in- 
cludes a summation over secondary hadrons in energy-energy 
correlations). In experiments, instead, the jet axis is determined 
exclusively from all particles inside the jet. The question of 



The effect of Varying t-mzn is more dramatic at MLLA. 
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In (k^/ 1QeV) 



n (kj^ / 1 QeV) 



Figure 5. CDF preliminary results for the inclusive distribution at 
various hardness Q in comparison to MLLA and NMLLA predictions 
at the limiting spectrum (Qo = ^qcd); the boxes are the systematic 



normalized to bin: ln(k^)=-0.1 (N') ! 

' CDF preliminary 




In (k^/ 1GeV) 

Figure 6. CDF preliminary results (Q = 119 GeV) for inclusive kj^ 
distribution compared with NMLLA predictions beyond the limiting 
spectrum. 



the matching of these two definitions goes beyond the scope 
of this letter. The NMLLA -spectrum has also been calcu- 
lated beyond the limiting spectrum, by plugging ( l24l i into ( |26] l. 
as illustrated in Fig. |6] However, the best description of CDF 
preliminary data is reached at the limiting spectrum, or at least 
for small values of A < 0.5, which is not too surprising since 
these inclusive measurements mostly involve pions. Identify- 
ing produced hadrons would offer the interesting possibility 
to check a dependence of the shape of k^^ -distributions on the 
hadron species, such as the one predicted in Fig.|6] Moreover, 
the softening of the -spectra with increasing hadron masses 
predicted in Fig. |6]is an observable worth to be measured, as 
this would provide an additional and independent check of the 
LPHD hypothesis beyond the limiting spectrum. This could 
only be achieved if the various species of hadrons inside a jet 
can be identified experimentally. Fortunately, it is likely to be 
the case at the LHC, where the ALICE 119} and CMS L20J 
experiments at the Large Hadron Collider have good identifi- 
cation capabilities at not too large transverse momenta. 



4 Conclusion 

To summarize, single inclusive fc^ -spectra inside a jet are de- 
termined including higher-order O (as) (i.e. NMLLA) correc- 
tions from the Taylor expansion of the MLLA evolution equa- 
tions and beyond the limiting spectrum, A 7^ 0. The agree- 
ment between NMLLA predictions and CDF preliminary data 
in pp collisions at the Tevatron is impressive, indicating very 
small overall non-perturbative corrections and giving further 
support to LPHD |3 |. The MLLA evolution equations for inclu- 
sive enough variables prove once more (see e.g. |5|) to include 
reUable information at higher orders than MLLA. 
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